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Abstract 



We demonstrate that decay angle correlations in r r+ and tt decay modes 
^ ■ could allow a determination of whether or not a neutral Higgs boson is a CP 
eigenstate. Sensitivity of the correlations is illustrated in the case of the e'^e^ 
\ ZH and H production processes for a two-doublet Higgs model with CP- 

violating neutral sector. A very useful technique for minimizing 'depolarization- 
factor' suppressions of the correlations in the ti mode is introduced. 



Determination of the CP nature of any neutral Higgs boson that is directly 
observed will be crucial to fully unravelling the nature of the Higgs sector. The 
! Standard Model (SM) Higgs boson is CP-even, while the Minimal Supersymmetric 
Model (MSSM) predicts two pure CP-even and one pure CP-odd neutral Higgs 
i. boson. More generally, it is entirely possible to have either explicit or spontaneous 

rS I CP violation in the neutral Higgs sector. Indeed, the simplest non-supersymmetric 

' two-Higgs-doublet model (2HDM) and the supersymmetric Higgs two-doublet plus 

singlet model both allow for Higgs mass eigenstates of impure CP nature.^ Here 
we shall focus on the 2HDM, in which CP violation results in three neutral states, 
ifj^i 2,3, of mixed CP character. 

Various types of experimental observables can be considered for determining 
the CP character of a given Higgs boson. The most direct probe is provided by 
comparing the Higgs boson production rate in collisions of two back-scattered- 

laser-beam photons of various different polarizations. A certain difference in rates 
for different photon helicity choices is non-zero only if CP violation is present, and 
has a good chance of being of measurable size for many 2HDM parameter choices. 
In the case of a CP-conserving Higgs sector, the dependence of the 77 
cross section on the relative orientation of the transverse polarizations of the two 
colliding photons may well allow a determination of whether a given H is CP- 
even or CP-odd!^'*' Note that the general utility of the photon-photon collision 
polarization asymmetries derives from the fact that the (one-loop) 77 couplings 



of a CP-odd and CP-even H are similar in magnitude, so that sensitivity is not 
strongly dependent upon the CP nature of the H. 

Correlations between decay products can also probe the CP nature of a Higgs 
boson. In this paper we focus on effects that arise entirely at trec-levcl, i.e. that 
do not rely on imaginary parts generated at one-loop. For the dominant two-body 
decays of a Higgs boson, we can define appropriate observables if we are able to 
determine the rest frame of the Higgs boson and if the secondary decays of the 
primary final state particles allow an analysis of their spin or helicity directions. 
An obvious example is to employ correlations between the decay planes of the decay 
products of WW or ZZ vector boson pairs and/or energy correlations among the 

decay products!^""' However, these will not be useful for a purely CP-odd H 

(which has zero trcc-lcvcl WW, ZZ coupling and thus decays primarily to FF) 
or for a mixed-CP H in the (most probable) case where the CP-even component 
accounts for essentially all of the WW, ZZ coupling strength (thereby yielding 
'apparently CP-even' correlations). In contrast, H decays to t~t+ or ti, followed 
by r or t decays, do, in principle, allow equal sensitivity to the CP-even and CP-odd 
components of a given Higgs boson. 

Indeed, a H eigenstate couples to FF according to: C ex F[a-\-ilry^)FH ,which 
yields 

{F+F+\H) o^h + ia(3F] {FJF_\H) (xh - ia(3F , (1) 



where (3p = 1 — Am^/rnjj, helicity-flip amplitudes being zero. The crucial point 

is that, in general, a and b are of comparable magnitude in a CP- violating 2HDM. 
In the notation of Ref. [12] we have 



Qjq '7-+ 



V sin P ' 

rrirCi 

V cos (3 ' 



K-T+ = - 



mt-siss cos l3 

rrirSiSs sin (3 
V cos P 



(2) 



where v — 246 GeV, tan /3 = t;2/i'i is the ratio of the vacuum expectation values of 
the neutral Higgs fields that couple to up and down-type quarks, and ci, si, C3, S3 
are cosines and sines of neutral Higgs sector mixing angles ai, a^; the couplings of 
Eq. (2) are those that appear in the Euler angle parameterization of the first Higgs 
eigenstate as defined in Ref. [12]. In this same notation, the WW coupling of H is 
ci cos f3 + S1C3 sin f3 times the WW coupling of the SM Higgs boson. Note that we 
have assumed a Type-II (as defined in Ref. [1]) 2HDM model, in which charged- 
lepton couplings of the H are analogous to down-quark couplings. It is consistent 
to take H to be the hghtest eigenstate and to allow the angles cti^s appearing in 
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ci,si,cs,ss to take on arbitrary values. Clearly, a and b are generally of similar 
size. SM couplings for the H are obtained in the case ai — (3 and as = 0. 

The use of azimuthal angular correlations in t~t'^ and tt decays to determine 
the CP eigenvalue of a pure CP state was explored in Ref. [4]; the additional 
azimuthal angle dependence that is present only for a mixed-CP cigcnstate has 
been noted in Ref. [8] in a special case of ti decays. Here, we shall present a unified 
treatment aimed at realistically evaluating the possibility of using correlations in 
H T~T^ and ti final states to determine if a decaying Higgs boson is a mixed 
CP eigenstate, thereby directly probing for the presence of CP violation in the 
Higgs sector. 

An efficient framework for our analysis is that developed in Refs. [13,4]. Con- 
sider the charged current decay F — > i?/, where F is a heavy fermion, / is a light 
fermion whose mass can be neglected, and R can be either a single particle or a mul- 
tiparticle state with known quantum numbers and, therefore, calculable coupling 
to the charged weak current. (Examples are r Ru, where R — tt, p, Ai, . . ., and 
t — >■ Wb, where W decays to a fermion plus anti-fermion.) For the i?'s of interest, 
the form of the hadronic current J^, deriving from the standard V — A interaction 
for the Jjj, = (i?|V^ — A^|0) coupling, is completely determined in terms of the final 
particle momenta. Using the particle symbol to denote also its momentum, and 
defining 

H^ = 4Re J^/ • r - 2f^J -r, 2e^p,,Im J^J*7^ (3) 

all useful correlations in H ^ FF decay can be obtained by employing the quan- 
tities 

a; = F • (H - H^), R^ = ml{U - H^)^ - F^F ■ (H - H^) , (4) 

and their F analogues. In the F rest frame, Rq — 0, R — m|.(n — H^), and 

|-R| = rriFtu. In fact, Sp = R/{mFUj) acts as an effective spin direction = 1) 

when in the F rest frame. 

Let us give some illustrative examples. For — 7r~z/ decay, oc vr^ and 

S*^- = 7r~ is the unit vector pointing in the direction of the 7r~'s three momentum 
(using angles defined in the t~ rest frame). For r~ p~z/ — 7r~7r'^z/, oc 

(tt" — 7r^)p, yielding H^ oc 4(7r~ — 7r*^)/x'^ ■ (ti"" — '^'^) + 2z/^mp, and, thence, Sp oc 

— * — * „ 

mr{7r~ — 7r'^){E^- — E^o) + where the pion energies and directions are 

defined in the t~ rest frame. For t — > W^b — > I'^vb, oc m(z/)7^(1 — ^^)v{l'^), 



■k For fixed tan (3 and w, there are 7 independent parameters in the 2HDM model in which the 
discrete symmetry that guarantees the absence of flavor-changing neutral currents is only 
softly-broken.'^^' These can be taken as (i = 1, 2, 3), the masses of the three neutral Higgs 
bosons, and the mass of the charged Higgs boson. 
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and oc l^u ■ b + z/^/+ ■ 6, 11^ oc • b — /^z/ • 6, so that 11^ — 11^ oc implying 

St = in the t rest frame. 

If the full (11 — n^)^ can be determined on an event-by-event basis, then we can 

define the 'effective spin' vectors Sp and for each event, and the distribution 
of the Higgs decay products takes the very general form 



dN(x 



{b^ + a^(5l){l + cos e cos 9) + (6^ - a^pij,) sin 9 sin 9 cos((/) - 0) 

-1 - - 

- 2a6/5i? sin 9 sin 6* sin(0 — (p) d cos 6*0? cos 9d(pd(p , 



where 6', and 6', define the angles of Sp and in the F and F rest frames, 
respectively, employing the direction of F in the H rest frame as the coordinate- 
system- defining z axis. (Note that the F coordinate axes are to be used in both 
the F and F rest frames to define the angles appearing in Eq. (5).) To employ 
Eq. (5) we must retain the ability to distinguish F from F. In fact, to determine 
sin(0 — 0) it is absolutely necessary that we be able to determine the F and F rest 
frames, i. e. their line of flight in the Higgs rest frame. We shall return to this issue 
momentarily. 

We note that at one-loop a and b can develop imaginary parts. In this case, 
a'^,b'^,ab should be replaced in Eq. (5) by |ap, |6p,Re(a6*), respectively. In addi- 
tion, new angular dependences arise in Eq. (5) of the form —2lm{ab*)j3F {cos 9 + 
cos 9). In principle, the imaginary parts of a and b are also sensitive to CP violation 
in the Higgs sector (but also to other types of CP violation — in the SM, non-zero 
effects appear at the 2-loop level). After including branching ratios, the statistical 
significance associated with isolating the above term is not very large if the 2HDM 

provides the only source of CP violation!"' In any case, these extra terms will not 
contribute to the correlations upon which we focus. 

If we cannot determine (H — H^)^ for each event, then Eq. (5) must be modi- 
fied. An extreme example is F — > Rf decay where the R decay products are not 
examined. In this case the angles of R in the F rest frame would be employed in 
Eq. (5), and 'depolarization' factors arise as a result of event averaging. In deriving 
Eq. (5), the angular independent term is actually multiplied by {mFUJF){rnFOJ-p) 

and the cos 9 cos 9, sva.9 sin ^ sin(0 — 0) and sin 9 sin 9 cos(0 — 0) terms by | 1 1 Rp] ■ 
On an event-by-event basis the ratio of these coefficients is unity, as outlined ear- 
lier. When averaged over events, this is no longer true. Consequently, when event 



-k Defining the 4- vector 5^ = Rfj,/{uimF), and similarly for 5^^, the underlying covariant form 

of the matrix clement squared is cx (a'^ + b'^){F -J^ - m%S -S) -\- {a^ - b'^){F - T^S ■ 

S -F ■ SF -S - Tn,%) - 2abeo,i:jpaS"^'^^ FPF" , with the convention eoi23 = +1- We note that 
this result reduces to that of Ref. [8] (except for a difference in the sign of the e term) in 
the case of tt where (using their notation I = Z+) = {mpli^/F ■ I — F^/niF)- 
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averaging (denoted by (...)) all the angle-dependent terms in Eq. (5) must be 
multiplied by Dp = (\Rp\) / {riiF{uF)) and/or its analogue, relative to the 
angle-independent term. We define D = DpDp. 

At first sight, the necessity of event averaging arises in the case of the ti final 
state, for which we will find that we must have one top decay leptonically and the 
other hadronically in order to define the ti line of flight and, thereby, appropriate 
angles in Eq. (5). For the hadronically decaying top, the problem is to distinguish 
the quark vs. anti- quark jet coming from the 1^ so as to construct (11 — H^)^ 
(which is proportional to the anti-quark (quark) momentum for t (i) 

decay) for each event.^ If we simply sum over all W decay product configurations, 
then the appropriate depolarization factor is easily computed by using oc 

and summing over W polarizations. One finds 11^ = and 11^ oc 6^ -|- W^{m'^ — 
m^r)/m'^. Employing this result yields Dt = {ml — 2m^)/(m^ + 2m^) ~ 0.4 
for mt = 174 GeV. In the modified Eq. (5) the angles for the one hadronically 
decaying t (or t) would then be those of the (or W~) in the t (or t) rest frame. 
(For the leptonically decaying t (or t) the angles of the l~ (or Z+) are directly 
measured and the associated Dj (or D^) is unity.) Similarly, if in r — > Ru the R is 
spin-1 and its decay products were simply integrated over, a depolarization factor 
of Dr = {m^ — 2m|j) / (m^ + 2m^) would enter. 

Fortunately, these severe depolarization factors can be avoided in both cases. 
For the bulk of r decays is a resonance of known quantum numbers decaying to 
easily distinguished particles, in which case we can construct (11 — 11^)^ event-by- 
event (see, e.g., the p example described earlier) and depolarization factors do not 
arise. In the case of a hadronically decaying t (t) , a simple helicity argument shows 
that the most energetic of the W'^ (W^ ) jsts in the t (i) rest frame is more likely to 
be the anti-quark (quark), i.e. the equivalent of the Z"*" (l^). Employing the angles 
of this most energetic jet (while integrating over the angles of all the other jets, so 
that the angles of this most energetic jet define the only direction associated with 
the decay) yields (via Monte Carlo calculation) D ~ 0.78, essentially independent 
of niH for niH ^ 800 GeV. 

Let us now specify our procedure for isolating the coefficients of the cos(0 — 0) 
and sin(0 — 0) angular correlation terms. Defining c = cos^, c = cos^, s = sinO, 
s = sin 6, Off) = cos 50, = sin 50 (where S(f) = (f) — (f)), and dJl = dcdcd5(f), and 
including a possible depolarization factor, we have 

I dN 1 . ^ _ _ _ . 

N~dn ^ + P2ssc^\ , (6) 



f We assume that there is no efficient technique for determining the sign of the charges of the 
quark jets resulting from the top decay. 
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where 



Pi = D 



+ a^(3l) ' 



2ahl3F 



D 



{h^ - a' Pi) 



(7) 



For a CP-conscrving Higgs sector, either a = or 6 = implying pi = and 
\p2\ = D. For a CP-mixed eigenstate, both a and h are non-zero. Thus pi 7^ 
provides an unequivocable signature for CP violation in the Higgs sector, while the 
difference D — \p2\ also provides a measure of Higgs sector CP violation. (Indeed, 
pi and p2 are not independent; Pi+ P2 — D^.) Values of pi ~ D and p2 ~ are 
common in an unconstrained 2HDM. 

To isolate pi and p2, we define projection functions /i^2(^; (^0) such that 
/ fi,2dQ = 0, / fi^2ccdfl = 0, J fisss^dfl = Stt, / fissc^dQ = 0, / f2SSS(j)dD, — 
0, and / f2SSC(j,dQ, — Stt. Then, pi^2 = / fi,2jj-'-^dfl. The critical question is 
with what accuracy can pi^2 be determined experimentally? In the absence of 
background, it is easily shown that the experimental errors of the determination 

arc given by 5pi,2 = (yi,2 - Pi,2)^''^/v^) where 1/1,2 = / /i,2l?S"^^' ^'^^ ^ 
total number of events. (For the /i,2 choices we shall make, yi^2 = / /i 2'^^/(87r).) 
If background is present, then this result is modified to 5pi,2 = [1/1,2 ^ Pi 2 + 
{B/S){yi^2 + P12 ~ 2pi,2pf2)]^'^^/A/5', where S is the number of H events, B is 
the number of background events, pf 2 is that for the background alone, and pi,2 

refers to the signal only. This result assumes that the background is precisely 
known, cither by detector Monte Carlo plus theory or high precision experimental 
measurement. The choices /i = (8/7r)e(s0) and /2 = (8/7r)e(c^) [where e(a) = +1 
(— 1) for a > (a < 0)] are equivalent to employing simple asymmetries, and yield 



(for B^O) yi,2 = (8/7r)2 and pi,2/5pi,2 = (7rpi,2/8) VA^/[1 - (7rpi,2/8)2] V2. (Note 



that (7r/8)pi,2 are the magnitudes of the (A^-^ — A'"_)/(A^+ + A^_)-type asymmetries.) 
For a functional form expressed in terms of orthogonal functions (upon integration 
over dQ), the error is minimized by using projection functions which match the 
angular dependence of the term of interest. Thus, we employ /i = (9/2)sss0 

and /2 = (9/2)ssc^, for which yi,2 = 9/2 and pi,2/(^Pi,2 = ^/2/9pl^2VN/[^ - 

(2/9)pf 2]^/^- Note that ^2/9 ~ 0.47 > (7r/8) ~ 0.39. 

We now discuss the Higgs production reactions and the r^r"*" and ti final state 
decay modes for which the angles of Eq. (5) can be experimentally determined. 
Consider first the t~t~^ case. The r decays are of two basic types: r luu 
and r — > Rv, where i? is a hadronically decaying resonance of known quantum 
numbers. Together these constitute about 95% of the r decays, with BR{t — > 

HRv) ~ 58.8%!^^' In the presence of two or more neutrinos, we cannot determine 
the r rest frame angles without employing a Higgs production reaction in which 
the Higgs rest frame can be determined without reference to its decays. Even 
if we know the Higgs rest frame, if either (or both) r's decay leptonically we 
will still not be able to determine either cos (50 or sin 50. Only by knowing the 
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Higgs rest frame and having both r's decay to Ru are there enough constraints 
to unambiguously determine For such decays it is crucial that the charge 

of R can be determined from an examination of its decay products; in fact we 
assume that the R decay products can be fully identified so that (11 — H^)^ can 
be determined for each event, thereby avoiding any depolarization factor. Thus, 
we employ D — 1 and an effective branching ratio for useful t~t^ final states of 
(0.588)^. Because of a two-fold ambiguity in the kinematic solution, determination 
of the sign of Scj) generally requires vertex tagging of both the r's!^^' Thus, the 
(0.588)^ should also be multiplied by the square of the efficiency for vertex tagging 
a r when estimating our ability to measure pi; however, because this efficiency is 
strongly detector-dependent we do not include it in our explicit numerical results 
for pi. 

In the case of ti decays, we cannot employ purely hadronic final states for which 
we would be unable to distinguish t from t. Even if the Higgs rest frame is known 
(and the four-momenta of both 6-jets are measured), double leptonic decays lead 
to a two-fold ambiguity in the determination of 50, so that only cos 50 could be 
computed. (As we see from Eq. (6), this is adequate for p2.) Only in the case where 
one top decays hadronically, and the other leptonically, are we simultaneously able 
to determine the exact ti decay axis and distinguish t from i. Thus, we employ an 
effective branching ratio for useful ti final states of 2 x (2/3) x (2/9) (keeping only 
I — e,ii). As noted earlier, employing one hadronic t (or i) decay and identifying 
the most energetic jet from the W'^ with the anti-quark (quark) leads to a 

depolarization factor of D ~ 0.78.^ Finally, we note that for a known Higgs mass 

and known top mass, the kinematical constraints in the single-leptonic ti final 
states are sufficient to determine unambiguously the momentum of the single i>, 
without knowing ahead of time the Higgs rest frame. This imphes that the ti final 
state correlations could in principle be employed at a hadron collider, although the 
extra initial and final state radiation present in hadronic collisions is very likely to 
lead to too much confusion for this to work in practice. 

In order to assess our ability to experimentally measure pi and p2, we have 
examined H production in the reactions e'^e~ ZH at a future linear e^e~ 
collider and p+p^ ^ iJ at a possible future p"'"p~ collider!^^' (In the ZH reaction, 
we employ both hadronic and I = e, n charged-leptonic decay modes for the Z, with 
net total branching ratio of ~ 76%.) For both reactions, the Higgs rest frame can 
be determined without reference to the H decays. For the e+e~ coUider we have 
adopted the optimal energy, ^/s — mz + V^mn, as a function of Higgs mass, and 



* The sin (5^ dependence involving just the two charged leptons in the double leptonic ti final 
state pointed out in Ref. 8 thus cannot be experimentally isolated. 

f If the Higgs rest frame is known, p2 can be obtained for double-leptonic decays with branch- 
ing ratio (2/9)^ and D = 1. Combining this channel with the semi-leptonic channel would 
result in roughly a 13% increase of the statistical significance values for a CP violation 
signal that we shall present in the case of p2. 
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assumed an integrated luminosity of 85 fb^ . Our results for statistical significances 
will assume that this ZH mode is essentially background free. However, we have 
not incorporated any efficiencies for cuts that might be required to guarantee this. 
We have also not included the dilution due to the ZZ continuum background for 
iTiH values in the vicinity of mz- Depending upon the detector resolution, this 
background can be STibstantial for mjj values between about 75 and 105 GeV. 
In this interval, our results should (at best) be considered an upper bound. For 
the [x^ [jT collider we have computed the Higgs signal and the continuum t~t'^ 
and ti backgrounds assuming unpolarized beams and a machine energy resolution 
of 0.1%, with ^/s centered at the (already known) value of itlh. We adopt an 
integrated luminosity of 20 fb^^. This is consistent with the hoped for integrated 

luminosity of 100 — 200 fb^^ of a multi-TeV collider when run at the lower 

energies required for direct production of the H in the mass range considered. 

For rriH values such that the e'^e~ — > ZH production mode is background 

free, the statistical significance of a non-zero result for pi is that given earlier, 

N^jj = \pi\/6pi, where Spi = (9/2 — pf)^/"^ /y/N, and is the number of events 
after including the branching ratios required to achieve the final state of interest: 
BReff = BR{H FF) x BR{FF X), where the latter FF branching ratios 
to useful final X states were specified above. In the case of p2 we must actually 
determine the statistical significance associated with a measurement of D — |p2|- 
This is given by iV|^ = [D - \p2\]/Sp2, where 6p2 = (9/2 - pl)^/^/^/N. 

In p^ p~ —>■ H, the continuum backgrounds must be included. The CP- 
conserving background does not contribute to pi, and the statistical significance 
of a non-zero value for pi is given by Ngj^ — \pi\/Spi with Spi — [9/2 — pi + 

(B/ S){9/2 + pI)]^^'^ /y/S, where S is the total number of events from H produc- 
tion, and B is the total number of events from the continuum background, in the 
final state of interest. Although the background does have substantial cos S(f) de- 
pendence, not only should wc have an excellent theory plus detector Monte Carlo 
simulation of the background, but also p2 for the background could be directly 
measured for ^/s values on either side of the Higgs resonance. Thus, we neglect 
errors in the background subtraction. In this case, we have N^jj = [D — \p2\]/5p2 

with 5p2 = [9/2 - pI + {B/S){9/2 + p^ - 2p2p^)]^/^/^, where pf is that for 
the background alone. (The non-superscripted pi 2 arc always those of the signal 
alone.) Results presented for Ngjj include p^, but differ negligibly from those 
obtained if is set to zero in the error expression above. 

Our results for the maximum N^jj and N^j^ values are presented in Fig. 1, 
where we have adopted a top quark mass of 174 GeV. The maximum values were 
found by searching over all values for the Euler angles ai and appearing in 
Eq. (2), holding tan/3 and mu fixed. In general, N^j-, is only slightly larger than 
Ngjy, as is easy to understand from the fact that \pi\ and D — \p2\ both have 
maximum values close to D. The production rates and branching ratios both 
depend upon the couplings of Eq. (2), as well as couplings to other fermions and to 
1^1^ and ZZ pairs. Couphngs to up-type fermions are, of course, analogous to the 
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Figure 1: The maximum statistical significances Ngj^ and Ngjj for H t ( ) and 

H ^tt { ), in e+e- ^ ZH {L = 85 fb"^) and /x+zu" H {L = 20 fb"^) production, 

after searching over all ai and values at fixed niH and tan/3. In each case, curves for the 
three tan/3 values of 0.5, 2, and 20 arc shown. In the t~t~^ {tt) mode Nsd values increase 
(decrease) with increasing tan/3, except in the case of H ^ tt, where the lowest curve 

is for tan/J = 0.5, the highest curve is for tan/3 = 2, and the middle curve is for tan/3 = 20. 

ti coupling given in Eq. (2) with mt replaced by the up-type fermion mass, while 
couplings to down-type fermions are analogous to the t~t'^ coupling of Eq. (2) 
with rrij- replaced by the down- type fermion mass. In computing the t~t^ and 
ti branching ratios, the full set of possible H decays to //, WW, and ZZ are 
included. As noted earlier, the results of Fig. 1 for Ngj^ in the t~t~^ mode do not 
incorporate efficiencies for r vertex tagging, required to determine the sign of S(f) 
(as needed for computing pi), and thus should be multiplied by the efficiency {not 
its square) for r tagging. Fortunately, it is expected that this efficiency will be 
relatively high for appropriate detector designs. 

Consider first the results for e~^e~ ZH collisions. From Fig. 1 we find that 
detection of CP violation through both pi and p2 is very likely to be possible for 
mji < 2mw via the H — >■ r^r"*" decay mode. This is an important result given that 
various theoretical prejudices suggest that the lightest Higgs boson is quite likely 
to be found in this mass range. For mu between 2mw and 2mf, a statistically 
significant measurement of CP violation will be difficult. For mn > '^rnt, detecting 
CP violation in the ti mode would require a somewhat larger L (of order 5 times 
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the assumed luminosity of L = 85 fb ^ for tan/3 between 2 and 5). 

In fi'^ — > H production, Fig. 1 shows that the maximum N^jj and A^^^^^ 

values in the t~t^ mode can remain large out to large Higgs masses if tan {3 is large, 
but that for small to moderate tan (5 values the statistical significances are better in 
e+e" ZH collisions when itlh < 2mw- The reason for this is obvious — Higgs 
production in fi^fi^ collisions is strongly enhanced at large tan f3. However, Fig. 1 
also indicates that the H ^ t^t^ channel has the advantage of possibly 

small sensitivity to the WW decay threshold at mu ~ 2mw- Such insensitivity 
arises when the Euler angles ai, as are chosen so as to minimize WW, ZZ couplings 
(and hence H — > WW, ZZ branching ratios) without sacrificing production rate. 

Thus, for L — 20 fb~^ A*"*"/^" coUisions could allow detection of CP violation all 

the way out to 2mt for tan (3 ^ 10. The ti final state extends the range of mn 
for which detection of CP violation might be possible only somewhat, and only 
if tan (3 lies in the moderate range near 2. A factor 10 higher n^iJ.~ luminosity 
{i.e. L = 200 fb~^, requiring a machine design focusing on center-of-mass energies 
below 1 TeV and, possibly, several years of running) would extend the range of 
possible detection in both modes: for tan/3 ^ 2, rriH values up to 2mt could be 

probed in the t~t~^ mode, while the ti mode might be useful out to quite high 
masses. 

Of course, in obtaining the above results we have implicitly assumed that 
the H does not have additional decays. If it is not the lightest Higgs eigenstate, 
decays of the H to a pair of hghter Higgs bosons or to Z plus a lighter Higgs 
boson might be kinematically allowed. If present, they would dilute the statistical 
significances of Fig. 1. However, if decays involving other Higgs eigenstates are 
significant, then there are many other direct 'signals' of CP violation in the Higgs 
sector that could be present. For example, simultaneous presence of H2 HiHi 
and either H2 ZHi decays (in our notation, H2 is the heavier state) or e~^e~ 
Z — > H1H2 production at a significant level would alone require the if's to be a 
mixture of CP-even (allowing decays to a pair of Higgs) and CP-odd (allowing Z 
plus Higgs decay) states. As another example, a substantial production rate for 
H12 in e~^e~ — > ZH12 combined with either the existence of H2 — * ZHi decays 
or e+e~ Z ^ H1H2 production would imply that the couplings ZZH\, ZZli2 
and ZH1H2 are all non-zero, which requires CP violation in the Higgs sector!^* 

We have not explicitly analyzed the case of a (non-minimal) supersymmet- 
ric model with CP violation in the Higgs sector. However, several comments are 
useful. First, decays to superpartner pairs (neutralino, chargino, slepton, squark 
pairs) might be important and would dilute the observability of pi and p2 in the 
T~T~^ and ti channels. In this case, one can consider using the the supersymmetric 



-k Strictly speaking, the above statements arc only true with regard to tree-level couplings; a 
ZZH coupling is present at one-loop even if the H is purely CP-odd. To completely avoid 
contamination from C-violating one-loop diagrams, three or more neutral Higgs bosons 
must be detected. For example, to all orders, non-zero values for all three of the couplings 
ZH1H2, ZH1H3 and ZH2H3 are only possible if CP violation is present. 
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partner pair events themselves. Generally, a measurement of pi,2 in superparti- 
cle pair channels would probe a subtle mixture of Higgs sector CP violation and 
CP violation deriving from complex phases in the soft-sup ersymmetry-br caking 
parameters that enter the chargino, neutralino, etc. mass matrices. However, re- 
strictions from neutron and electron electric dipole moments suggest that the latter 
phases are quite small, in which case an observable non-zero value for pi or D— \p2\ 
would imply CP-violation in the Higgs sector. Procedurally, those events in which 
each member of the supcrparticle pair decays to jets and/or charged leptons plus 
a single lightest neutralino {i.e. the LSP, of presumably known mass) would allow 
treatment along the same lines as the t~t~^ mode after correcting for the finite 
mass of the single invisible LSP. Generally, lifetimes would be too short for vertex 
tagging, and only \5(f)\ could be determined, so that only measurement of p2 would 
be possible. Of course, if the decays of a supersymmetric model H were spread 
out over many channels, there might be inadequate statistics in any one channel. 

It is also useful to comment on how well pi and p2 can be measured in the 
case of a CP-conserving Higgs sector. Recall that pi = and p2 — —D for a 
CP-odd, CP-cvcn Higgs boson. As an example, consider e~^e~ — > ZH, where H is 
CP-even. (Only a CP-even Higgs boson will have usable e~^e~ ZH production 
rate.) For simplicity, we assume that the H has SM couplings. (The results will 
then also apply to the Hghtest CP-even MSSM Higgs boson in the case where 
the other Higgs bosons have masses m^o ~ mjjo ~ TnH± ^ 2m^, in which limit 

the is SM-like.) For L = 85 fb~\ we find that the r^r+ mode yields 5 pi ~ 5p2 
increasing from ~ 0.05 to ~ 0.13 as mn ranges from 60 GeV up to 160 GeV. 
Comparing to 1 {D — I'va. this case), we see that a simultaneous measurement of 
pi and p2 would provide a very strong confirmation that a SM-like Higgs boson 
is indeed CP-even, pi would be particularly valuable for mu ^ Taz since the ZZ 
continuum background yields non-trivial cos 50 dependence. For mu > 2mw, even 
a rough determination of pi^2 would not be possible in this mode due to the rapid 
fall in the H t~t~^ branching ratio resulting from the onset of WW and ZZ 
decays. In the ti mode, dpi ~ dp2 ranges between 0.4 and 0.8 for rriH between 2mt 
and 750 GcV. Comparing to D ~ 0.78 we see that statistics would be inadequate 
to clearly distinguish between a SM-like CP-even Higgs boson and one of mixed- 
CP nature on the basis of the azimuthal angle correlations. Of course, in a general 
2HDM, the H in question might be CP-even but have reduced WW, ZZ couphng. 
(In this case, the reduced production rate via e~^e~ ZH would be apparent, 
but would not on its own indicate whether the Higgs was CP-even or of mixed-CP 
character.) For reduced WW, ZZ coupling, the H WW, ZZ branching fractions 
decline much more rapidly than the ZH cross section and the errors on pi and p2 
in the tt mode can be sufficiently small to provide a strong indication of the CP 
character of the H. 

In summary, our results show that if the Higgs sector is CP-violating then 
there is a substantial possibility of explicitly exposing this CP violation through 
azimuthal angle correlations between final state particles m H ^ t^t^ , where the 
H is produced via e+e~ ZH or p+p^ — * H, with assumed integrated luminosi- 
ties of L = 85 fb~^ and L = 20 fb~^, respectively. Of particular importance is the 
general utihty of the r^r"*" mode in e+e" — > ZH coUisions for Higgs masses in the 
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theoretically preferred mn ^ 2mw region. For mn > '^^mt, azimuthal correlations 
in the ti mode also provide sensitivity to CP violation. However, statistically reh- 
able correlation measurements arc predicted to be possible for a smaller portion of 
parameter space, which, however, would expand considerably if higher luminosi- 
ties were available. The correlations employed rely on the fact that CP violation 
generally leads to non-trivial dependence (not present if the Higgs sector is CP- 
conserving) on the sine of an appropriately defined azimuthal angle (50), and to 
dependence on the cosine of Scf) that is substantially different than that predicted 
when CP is not violated (see Eqs. (5) and (6)). Of the two possible CP-violation- 
sensitive correlations, pi (obtained from the sin 50 dependence) provides the best 
opportunity for detecting CP violation in the 2HDM Higgs sector, both because of 
a somewhat larger statistical significance, and because the associated non-trivial 
dependence on sin S(f) cannot arise from CP-conserving backgrounds or detector ef- 
ficiency effects. However, p2 (deriving from dependence on cos 50) provides nearly 
as good a probe of CP violation. Further, there is a tendency for both CP violation 
'signals' (namely a large value for pi and a small value for P2) to be simultaneously 
substantial as a function of the two-Higgs-doublet model Higgs mixing angle pa- 
rameters. On the experimental side, measurement of pi in the t~t^ mode requires 
high efficiency for r vertex tagging, in order to determine the sign of the crucial 
50 azimuthal angle, while prospects for measuring both pi and p2 in the ti mode 
could be improved somewhat if a still more efficient technique for identifying the 
quark (anti-quark) jet in W~ {W~^) decay were available, e.g. by determining the 
charges of the jets coming from a hadronically decaying W. 
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